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In order to elucidate the vacuum structure of the Aoki phase, we carried out a numerical investiga- 
tion of QCD with two flavours of Wilson fermions, within the p.d.f. framework and in the absence 
of external sources. The simulations performed at V = 4 4 suggest a rich vacuum structure, where 
the observable ixfrys V is allowed to take non-zero values of the same order of magnitude than the 
order parameter of the Aoki phase ixjfys^Y- However, the simulation at higher volumes V = 6 4 
suffers from large statistical errors. 
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1. Introduction 

A year ago, we applied the p.d.f. formalism [|l|] to QCD with 2 flavours of Wilson fermions, 
and in particular to the Aoki phase ph, to analyze its spontaneous flavour and Parity breaking pat- 
tern [Q]. The investigation led us to several interesting and undocumented results, which depended 
on the properties of the spectral density pu (ju) of the Hermitian Dirac Wilson operator H = 75 A. 
The spectral density should be an even function of the eigenvalues in the continuum limit a — > 0, 
but it has been also suggested that it may become symmetric already in the thermodynamic limit 
V — > °°. These two different behaviours come out into different possibilities: 

1. If the spectral density becomes an even function of the eigenvalues only in the continuum 
limit, then violations of hermiticity are expected to induce large artifacts in the measurement 
of the T] —meson mass at finite lattice spacing, even if V — > °°. 

2. If, on the contrary, it is enough to reach the thermodynamic limit to find a symmetric spectral 
density, that is the most plausible case [J|, Q], then two new scenarios appear as V — > °°, where 
the Aoki phase displays a rich, unexpected behaviour. In the first scenario, the conditions 

(n/0w) 2 ")/0 «GN (1.1) 

are verified as a consequence of the spectral symmetry [|p. The standard picture of the 
Aoki phase is incomplete, and new phases, unrelated to the original Aoki phase, appear. 
The second possibility, predicted by the chiral lagrangians |J], complies with the standard 
picture of the Aoki phase. The following equations 

(iy 75¥ ) 2 ")=0 «£N (1.2) 



are satisfied, and this give rise to an infinite set of independent sum rules for the eigenvalues 
/l, one for each n. These scenarios are mutually exclusive. 

Up to now, there is no theoretical proof selecting one of these possibilities. In [Q], an antisymmetric 
spectral density in the thermodynamic limit was excluded, based on quenched simulations. There- 
fore we decided to perform the dynamical fermions simulations required to distinguish between 



(jy) and ( ]L2p . 

The current paper is organized as follows: In the next section, the details and technical dif- 
ficulties of the simulation of dynamical fermions inside the Aoki phase are reviewed. Section 3 
shows and analyses our numerical results, and the last section is devoted to our final conclusions. 

2. Simulating the Aoki phase without external sources 

Although the simulations of the Aoki phase with dynamical fermions are nothing new in the 
lattice QCD panorama, these have always been performed under very special conditions: an exter- 
nal source is added to the action with a twisted mass term 
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in order to (i.) regularise the small eigenvalues of the Dirac Wilson operator, which appeal - only 
in the Aoki phase and usually spoil any attempt of simulation without the external source, and 
(ii.) to analyse the pattern of spontaneous flavour and Parity breaking. As proved in the use 
of an external source like ( |Q| ) selects a standard Aoki vacuum, complying with (|Q|). The only 
way to investigate the existence of the new phases characterized by ( |0| ) is to remove the external 
source, and extract results from direct measurements of the Gibbs state (£-regime). The latter point 
is solved by the p.d.f. formalism, but the former -the removal of the external source- has been 
an unexplored option in the Aoki phase dynamical simulations for a number of years. The reason 
is the appearance of small eigenvalues of the Hermitian Dirac Wilson operator, of order 0[y). 
As far as we know, the technical problems inside the Aoki phase are similar to those faced when 
trying to reach the physical point: The critical slowing down spoils the efficiency of simulations, 
the Dirac Wilson operator becomes increasingly harder to invert, and the performance decreases 
dramatically. In fact, for some values of the parameters (j8 , fc) in the coupling-mass phase diagram, 
the standard algorithms to invert the Dirac operator just will not work. This fact called for a research 
on competent algorithms to simulate inside the Aoki phase. 

Inspiration came in the recent results for new algorithms which reduce the critical slowing 
down for small masses in several orders of magnitude. In this work, we successfully implemented 
a SAP preconditioner in a GCR inverter, as explained in The new inverter allow us to perform 
simulations inside the Aoki phase without external sources at a reasonable speed. Unfortunately 
this is not the whole history. 

The simulations done outside the Aoki phase feature a spectral gap around the origin, with 
a symmetric spectrum 1 . This gap is required to preserve Parity and flavour symmetry How- 
ever, the Aoki phase breaks both of them, so this gap was not present in our simulations, and as 
explained above, the smallest eigenvalues of each configuration took values of order 0[y). As the 
eigenvalues approach the origin, it may happen that they try to change sign, rendering the spectrum 
asymmetric, but this movement is forbidden by the Hybrid Montecarlo dynamics, i.e., the HMC 
algorithm is not ergodic for Wilson fermions inside the Aoki phase. Therefore we are introducing 
artificial constraints in the dynamics, and the final results are bound to be modified. That is why we 
considered another dynamical fermion simulation algorithm, the Microcanonical Fermionic Av- 
erage (MFA) algorithm |^|, which solves the problem of the eigenvalue crossing, but converges 
poorly as the volume increases. 

The problem was overcome by using an argument developed in ^ and applied to the current 
case in [Q], where the asymmetry of the spectrum n^ sym is related to the topological charge Q as 

nAsym x Q- (2-2) 

Then the constraints imposed by the hybrid Montecarlo are equivalent to leaving the topological 
charge fixed. Since the measurement of the observables should not depend on the value of the 
topological charge in the thermodynamic limit, we can select the symmetric sector Q = and 
measure our observables there, where one expects to have smaller finite volume effects. 
Another interesting possibility is (i.) to measure the weight of the different TiAsym 

sectors in the 

partition function via the MFA algorithm, (ii.) then perform Hybrid Montecarlo simulations within 
'By symmetric here we mean that the number of positive and negative eigenvalues are the same. 
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the relevant sectors, and (iii.) do a weighted average of the observables obtained in the HMC using 
the MFA weights. The results obtained with this method were contrasted to those coming from the 
direct MFA simulations, and with the HMC simulations at fixed nA sym . 



A third proposal for simulations was announced during the talk: The addition of a small (2.1) 
external source, large enough to regularise the small eigenvalues, but small enough to avoid vacuum 
selection, would allow us to include the eigenvalue crossing phenomenon in our HMC algoritm. 
Unfortunately, there exists no value of the external field h capable of this two deeds at the same 
time. If the field is very small, no vacuum is selected, but the eigenvalues do not cross the origin 
during the simulations; on the contrary, for larger values of the field h, the standard Aoki vacuum 
is selected. This third possibility was, thus, forsaken. 



3. Numerical results 

We performed measurements of three observables of interest 

Table 1: Expected behaviour of the analysed observables in the different scenarios as V — > °°. 





Outside 
Aoki 


Aoki 
Standard Wisdom 


Aoki 
Our Proposal 


({iVuYsVu?) 


~o 


^0 








~0 




((^75W) 2 ) 


~0 


7^0 





Our measurements refer always to the second moment of the p.d.f., which should be zero 
in case of symmetry conservation, and non-zero if the symmetry is spontaneously broken [jjj]. 
The first observable signals Parity breaking, and should be non-zero inside the Aoki phase. The 
second one allows us to distinguish between our proposal -there is an additional Aoki-like phase 
verifying ^(/i//y 5 i//) 2 ^ ^ 0- and the one involving an infinite number of sum rules. Finally, the third 
observable is the landmark of the Aoki phase, marking spontaneous flavour and Parity breaking. 

The first set of simulations were performed using the HMC algorithm, improved with a SAP- 
preconditioned solver. The symmetric runs riAsym = were performed starting from a cold (ordered, 
links close to the identity) and from a hot (disordered, close to strong coupling) configuration, 
obtaining identical results within errors. We could not find an asymmetric state riAsym = 1 in a 
cold configuration at the values of ?c explored, the asymmetric run was started only from a hot 
configuration. 

We expect all the observables to have non-zero expectation values, even outside the Aoki 
phase, due to finite volume effects. However, the values inside the Aoki phase are an order of mag- 
nitude larger than those measured outside the Aoki phase. Outside the Aoki phase, the following 
approximate rule holds: 

2((W5V«) 2 ) « ((iYYsY) 2 ) ~ ((iYYs^Y) 2 ) , 
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Table 2: Results of the Hybrid Montecarlo measurements. 





HMC V = 4 4 
Outside Aoki* 


HMC** V = 4 4 

f^Asym — 


HMC V = 6 4 

f^Asym — 


HMC V = 4 4 

H-Asym — 1 


NConf 


20002 


19673 


664 


10002 




2.098(3) x 10~ 3 
4.149(4) x 10~ 3 
4.244(4) x 10~ 3 


1.90(3) x 10~ 2 
2.69(12) x 10~ 2 
4.90(9) x 10~ 2 


9.2(30) x 10~ 3 
-4.3(30) x 10- 1 
4.7(30) x 10- 1 


6.52(72) x 10~ 3 
-4.49(50) x 10~ 2 
7.10(30) x 10~ 2 



*Point outside the Aoki phase j5 = 3.0, K = 0.22. 
**Point inside the Aoki phase j3 = 2.0, K = 0.25. 



which is a manifestation of (i.) the presence of an spectral gap and (ii.) the high level of symmetry 
of the spectral density, even at small volumes. This facts can be seen in the p.d.f. expressions for 
these observables in terms of the eigenvalues jx 



(i¥uY5YuY 



(iYYsY) ) = 



(iYYs^YY 




2 N 



(3.1) 

(3.2) 
(3.3) 



where the term (jjj -^j almost vanishes outside the Aoki phase. Nonetheless, inside the Aoki 
phase the gap disappears, the asymmetry near to the origin can become relevant and this rule would 
break down. 

The numerical results are reported in Table 2. Some comments are in order: inside the Aoki 
phase, due to the parity and flavour breaking, we expect to find (iYuYsYu) 2 an d (iYY5^¥) 2 to be 
large with respect to the case of standard QCD (outside the Aoki phase), and this is exactely what 
we get. The crucial point to discriminate between the two scenarios depicted above is the behaviour 
of (i\j/j5Y) '■ if we ta k e the results of nA Sym = in the 4 lattice at face value we are induced 
to conclude that the non-standard scenario for the Aoki phase is favoured, being the (Hj/ysy) 
expectation value of the same order of magnitude of the other two observables (and an order of 
magnitude larger than outside the Aoki phase). The results of the 6 4 lattice seem to add no useful 
informations due to the large statistical errors (we hope to get better quality results in the future). 
On the other hand the result for n^ym = 1 can seem strange at a first sight (a negative number for 
the expectation value of the square of an hermitian operator), but we have to take into account that 
we are restricting ourselves to a single topological sector and what is relevant is the relative weight 
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of the various topological sectors to the final result (that should be positive). The evaluation of the 
weights can not be performed using HMC simulations due to ergodicity problems. 

We also observed that, in order to achieve in the asymmetric run the same acceptance ratios 
(~ 90%) than in the symmeUic runs, we had to reduce the simulation step by a factor of ten, for 
the forces inside the HMC became much larger than expected. We took this as an indication that 
the system was trying to return to the symmetric state, pushing the eigenvalues through the origin, 
thus increasing in an uncontrolled manner the norm of the inverse of the Dirac Wilson operator. In 
the symmetric run, the eigenvalues certainly tried to cross the origin, however this did not happen 
continuously, but only from time to time. Hence, although we can not fully rely in any of these 
results because of the aforementioned problems of the HMC, we find the results of n^sym 

= much 

more believable of those of the asymmetric state. 

Then we introduced the MFA algorithm^]], with the hope of solving the ergodicity problems. 
As in the MFA the contributions of fermions is added during the measurement, the eigenvalues can 
cross the origin at will. So the MFA algorithm allowed us to measure the weights of the different 
sectors n^ sym = 0,1..., and then use these weights to correctly average the HMC data. 

Table 3: Weights of the different sectors according to MFA algorithm. 



Volume 


fl-Asym — 


ftAsym — 1 


fl-Asym — 2 


4 4 


85.9 ±5.8% 


14.1 ±5.8% 


0% 



Table 4: Results for weighted HMC. 





Weighted HMC 
V = 4 4 


{(iYuYsVuf) 
({iYYsY) 2 ) 
(WsW) 2 ) 


1.72(9) x 10~ 2 
1.7(6) x 10~ 2 
5.2(3) x 10~ 2 



4. Conclusions 

We have performed dynamical simulations of QCD with two flavours of Wilson fermions 
inside the Aoki phase, without any external source. We succeed overcoming the critical slowing- 
down of the simulations by using new algorithms developed recently, but our data showed that the 
HMC algorithm is not ergodic inside the Aoki phase without twisted mass term. The measurements 
collected reveal the presence of a Parity and flavour breaking phase through the non-zero expecta- 
tion value of the operators (i^uYsWu) 2 and (iYuY5^3Wu) > where the eigenvalues of the Dirac Wil- 
son operator could take small values up to O (A). Concerning the expectation value of (/y/ysy/) 2 , 
the results in the 4 4 lattice point clearly in favour of a non-standard scenario for the Aoki phase, 
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whereas the 6 4 lattice results are not conclusive because of larger statistical errors. The possibility 
of increasing the statistics and the volume of the lattices used in the simulations would require large 
computer facilities, due to the high cost of the inversion of the Dirac Wilson operator in presence 
of small eigenvalues. Therefore we are working in theoretical arguments that may help to clarify 
the issue. 
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A. Notation 

Even though the p.d.f. is quite a powerful tool, its use is not very widespread in quantum 
field theories. That is why the language of the p.d.f. is a bit different to that of QFT's. This fact 
sometimes leads to confusion and misunderstandings, therefore we want to devote a few lines to 
explain the notation used thorough this paper. 

When we are referring to the observable iyjsy, what we mean is 
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Thence, the expectation value of the operator is computed as 




and its higher powers become 
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